First-principles theoretical evaluation of crystalline zirconia and hafnia as gate oxides 

for Si microelectronics 



(N 
O 

o 

(N 

< 



Vincenzo Fiorentini and Gianluca Gulleri 
INFM and Dipartimento di Fisica, Universita di Cagliari, 
Cittadella Universitaria, 1-09042 Monserrato (CA), Italy 
(Dated: February 1, 2008) 

Parameters determining the performance of the crystaUine oxides zirconia (Zr02) and hafnia 
(Hf02) as gate insulators in nanometric Si electronics are estimated via ab initio calculations of the 
energetics, dielectric properties, and band alignment of bulk and thin-film oxides on Si (001). With 
their large dielectric constants, stable and low-formation-energy interfaces, large valence offsets, and 
reasonable (though not optimal) conduction offsets (electron injection barriers) , zirconia and hafnia 
appear to have a considerable potential as gate oxides for Si electronics. 

PACS numbers: 68.35.-p, 77.22.-d, 85.30.-z, 61.66.-p 
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The performance needs of modern information tech- 
nology are forcing Si-based ultra-large-scale-integrated 
(ULSI) devices into the domain of nanometric dimen- 
sions. This downscaling implies, among others, the ef- 
fective continuing reduction of the physical thickness of 
insulating gate oxide layers in CMOS (Complementary 
Metal-Oxide-Semiconductor) devices. Amorphous Si02, 
the natural oxide of Si technology, is now nearing its fun- 
damental size limits, with physical thicknesses currently 
down to 2 unit cells This leads to uncomfortably 
large (> 1 A/cm^) leakage currents and increased fail- 
ure probabilities. The main reason for the strong reduc- 
tion of gate-oxide thickness in device downscaling is the 
need for increasing capacitances in the CMOS conducting 
channel. In a CMOS, the gate oxide layer dominates the 
series capacitance of the channel. An increase in capac- 
itance can be obtained reducing the dielectric thickness 
d/e of the oxide layer, having physical thickness d and 
relative dielectric constant e. Given its small dielectric 
constant, it is understandable that Si02 as a gate ox- 
ide has emerged as one of the key bottlenecks in device 
donwscaling |2| . 

It thus appears that, if Moore's law Q on ULSI circuit 
component density - and hence circuit performance - is 
to remain valid in the next decade, a replacement will 
have to be found for silica as a gate insulator. The ba- 
sic selection criteria for such a replacement are i) larger 
dielectric constant ("high-K"), ii) interface band offsets 
to Si as large as or comparable to those of silica (es- 
pecially the electron injection barrier), in) epitaxy on 
Si energetically not too costly, iv) thermodynamical sta- 
bility in contact with Si. In this work we address the 
expected performance, in terms of the above criteria, 
for the two important current candidates ||l|, |[ ^ haf- 
nia (Hf02) and zirconia (Zr02) through first-principles 
density-functional calculations of the structure, energet- 
ics, thermodynamical stability, dielectric constants, and 
band offsets of crystalline hafnia and zirconia thin films 
epitaxially grown on the (001) face of crystalline Si. We 
find stable and moderate-cost interfaces, large dielectric 



constants, and large band offsets, except for the electron 
injection barrier, estimated at 1 eV at most, appreciably 
lower than the Si/silica barrier. 

Our density functional theory calculations in the gen- 
eralized gradient approximation |^ use the VASP code 
1^ and the ultrasoft pseudopotentials provided there- 
with. Semicore states are treated as core for Hf and Zr; 
test calculations done including the semicore as valence 
using the all-electron PAW Q method as implemented 
in VASP joj confirmed the pseudopotential results. Bulk 
optimizations were done in a 12-atom (conventional fee 
or fct) cell, while the interfaces are simulated by (001)- 
oriented oxide/Si superlattices contained in tetragonal 
cells of c(2 X 2) basal section, and in-plane lattice con- 
stant asi=5.461 A, our theoretical value for bulk Si. In- 
terface supercells contain around 50 atoms depending on 
the local interface structure, with 9 layers (18 atoms) for 
the Si region, and typically 11 layers (e.g. 24 oxygen and 
10 Zr atoms) for the oxide region. The plane-wave ba- 
sis cutoff is 350 eV; for the k-space summation we use 
4x4x4 meshes for the bulk and 4x4x1 meshes for the 
z-elongated interface supercells. 

Bulk and Si-epitaxial structure - Bulk hafnia and 
zirconia were studied in the fluorite, monoclinic, and Si- 
epitaxial structures. The lattice parameters for Zr02 
are a— 5.10 A for fluorite, and (a,&,c) = (5. 186,5. 255,5.351) 
A, off-normal angle 0=8.83° for monoclinic. For Hf02, 
a=5.06 A for fluorite, and (a,6,c) = (5. 108,5. 175,5. 280) 
A, off-normal angle 0=8.80° for monoclinic. The latter 
phase is favored over fluorite by 0.115 eV/formula unit 
for Zr02 and by 0.248 eV/formula unit for Hf02. The re- 
sults agree with experiment and with recent calculations 
The formation enthalpies AHq 
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are -11.52 

eV and -10.74 for hafnia and zirconia respectively (close 
to experiment, as usual using GGA) compared to -8.30 
eV for silica: therefore both oxides are stable in contact 
with Si with respect to the decomposition into silica and 
metal. The same holds for the epitaxial phase discussed 
next, whose excess energy is only about 0.2 eV/formula 
above the monoclinic. 



2 



The tetragonal Si-epitaxial crystalline phase of each 
oxide was obtained imposing the in-plane lattice constant 
of Si, and adjusting the axial ratio and internal coordi- 
nates in the 12-atom conventional cells. The axial ratios 
c/asi are 0.92 for Zr02 and 0.90 for Hf02. We verified 
by variable-cell damped dynamics Q that this tetrago- 
nal bulk is stable against monoclinic distortions. The Si- 
epitaxial configuration, depicted in Fig. |l]for Zr02, may 
be viewed as a z-stacking of cation-anion bilayers alter- 
natingly oriented at 90° to each other, in which a) metal 
cations are disposed in dimerized (llO)-like rows (cation- 
cation distances within the rows 3.4 and 4.2 A compared 
to 3.86 A ideally), and b) oxygens quadruplets, originally 
square in fluorite, elongate to rhomboids along the (110) 
rows bending slightly sideways. The cation (anion) co- 
ordination decreases from 8 to 6 (from 4 to 3), in partial 
analogy to the monoclinic structure |]l0| . 




FIG. 1: Si-epitaxial structure of Zr02. Grey (black) atoms: 
O (Zr). 

The elastic energy i^^jastic Si-epi distorted bulk 

is 0.23 eV/formula or 5.87 meV/A^ for Zr02, and 0.16 
eV/formula or 4.37 meV/A"^ for Hf02 with respect to 
monoclinic bulk (i.e. both are slightly favored energeti- 
cally over fluorite, whose occurrence is anyway barred by 
symmetry). While substantial, these energies are com- 
parable to those of order ~4 meV/A'^ involved (for much 
smaller strains) in nitride semiconductor epitaxy |]l^ . As 
we now discuss, the knowledge of the volume-specific epi- 
taxial strain energy enables us to extract an area-specific 
interface energy, as well as to estimate the critical pseu- 
domorphic growth thickness. 

Interface energetics and offsets - Assuming a c(2x2) 
basal section, we investigated for both materials several 
local structures and terminations of oxide/Si (001) in- 
terfaces, e.g. Si/0, Si/metal, Si/metal-bilayer, mixed 
Si-metal layer/0, mixed Si-metal layer/0 with 50% va- 
cancies. The starting configuration of the oxide portion 
of the interface superlattices is assembled using the opti- 
mized Si-epi structure. The supercell length and atomic 



TABLE L Formation energies (eV/A^) of, and valence and 
conduction band offsets (eV) at different Si (001)/oxide in- 
terfaces. The assumed growth conditions are indicated. The 
best offset /energetics combinations are displayed in under- 
lined bold for metal-rich conditions, and bold for oxygen rich 
conditions. AH GW corrections are included. 



Material 


— > 


HfOa 


Zr02 


Interface | 


Growth 


VBO GEO Efoon 


VBO CBO Efoon 


Si/0 
Si-M/O 
Si/M 
Si-M/O vac 
Si/0 vac 


O-rich 
stoich 
M-rich 
M-rich 
stoich 


4.14 0.47 -0.16 
4.40 0.19 0.17 
3.96 0.65 0.12 

4.62 -0.01 0.22 


4.08 0.72 -0.21 
4.18 0.62 0.12 
4.72 0.08 0.07 
3.91 0.89 -0.15 
3.70 1.10 0.13 



positions are then reoptimized: the axial ratio remains 
unchanged, and relaxations occur only in the first two 
interface-neigboring layers. The interface energy can be 
expressed as the difference of the energies i?sL of the 
interface cell, and £^buik of the corresponding bulk com- 
ponents, as 

Eform = ^[EsL - £^bulk] ^ ^[{2 A S + nsi^^Si^^Si + 
+ "-oxVox^^ox) - (»^SiVsi-Esi -I- UoxV^^E^^)] S 

with n the number of bulk units, V, V and E, E' the 
corresponding volumes and energies per unit volume, A 
the basal superlattice area. The formation energy per 
unit area, 5, can be extracted unambiguously if the oxide 
bulk energy is calculated in the same strain state as in 
the superlattice (Si remains unstrained), as in that case 
all volume-dependent terms drop off. Any other choice 
of the bulk energies inserts a volume dependence in the 
interface energy p^ . 

The interface cell may be stoichiometric, metal- or 
oxygen-deficient depending on its local structure. Its 
formation energy will therefore depend on growth con- 
ditions, metal-rich ones favoring oxygen deficit, and O- 
rich favoring oxygen excess. Theoretically, this is de- 
scribed by fixing the chemical potentials of the con- 
stituents. Here, only one potential - e.g. oxygen's - 
is independent: — ^02/2 nieans O-rich conditions, 
and fiQ = fJ.02/2 + AHox/2 metal-rich ones. 

The formation energies of the various interfaces are 
listed in Table |. The standard Si-0 interface is favored 
in O-rich growth conditions. In metal-rich conditions, 
the preferred structure is the mixed Si-metal to 50 % 
vacant oxygen layer interface depicted in Fig. ^, which 
remarkably is the same as was recently obtained [T^ in 
all-electron ab-initio molecular dynamics simulations of 
metal deposition on, and oxidation of. Si (001). Notably, 
the two favored interfaces have large negative formation 
energies (referred, we remind, to the pre-strained bulk). 
This energetic gain in interface formation will be coun- 
terbalanced by the excess energy of the film's upper sur- 
face, and by the build-up of epitaxial elastic energy in 
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the growing layer. An estimate of the critical thickness 
tc for pseudomorphic growth over an area A then results 
from 

AEform + Ate ^elastic + ^^^-rf = 0, 

using which we predict that crystalline zirconia and haf- 
nia thin films should grow pseudomorphically on Si (001): 
indeed, using our calculated values for, e.g., zirconia, 
and the GGA surface-energy estimate for the tetrago- 
nal phase Esurf ^ 0.05 eV/A^ [|l||, we obtain --18 A 
and 27 A for metal- and oxygen-rich conditions respec- 
tively. The poly-Si gate, forming a Si/oxide interface in 
the place of a free oxide surface, should further stabilize 
the structure. 




FIG. 2: The [mixed metal-Si]/[0 50% vacant] interface (black: 
Si; grey: metal; white-dotted: O) . 

The interface band offsets are evaluated for each inter- 
face using the standard 'bulk-plus-lineup' procedure [ p5[ , 
expressing the valence offset (VBO) as the sum of the in- 
terface potential lineup and the valence-band-top differ- 
ences of the separately-considered bulks. The conduction 
band offsets, hence the electron injection barriers, is es- 
timated as CBO = F^ll^" - E|ip - VBO. The gap of Si 
is taken to be 1.1 eV; for both oxides, we use our GGA 
gaps corrected with the GW data of Ref. (l^ for Zr02, 
namely 5.9 eV and 5.7 eV for zirconia and hafnia respec- 
tively. These values are close to experiment for hafnia, 
and near the bottom of the (large) experimental range for 
zirconia. We neglect spin-orbit corrections, which should 
be well below ±0.1 eV as the valence states are oxygen- 
derived. We do include, instead, the quasiparticle cor- 
rections to the bulk valence-band edges at the GW level: 
this is essential since these corrections are of order '-^1 
eV in oxides compared to typical ~0.1 eV in semicon- 
ductors. We apply to the VBOs an overall correction 
of -1.08 eV, resulting from the -0.15 eV correction ||T^ 
for Si and the -1.23 eV correction Q for Zr02. Using 
the latter for both oxides introduces some uncertainty 
in the Hf02 results, but unfortunately no GW data are 
currently available for hafnia. 

In Table | we report the predicted VBOs and CBOs. 
Qualitatively, VBOs cluster around 4 eV, with apprecia- 
ble structure dependence, and CBOs are in the range 



to 1 eV. Interestingly, for zirconia the energy-wise most 
favorable structures have some of the largest conduction 
offsets. The high-end CBOs, ~ 1 eV, are smaller than, 
but comparable to, the 1.4-1.5 eV estimates by Robert- 
son ijl^ , who used a simple charge-neutrality-level model 
at the empirical tight-binding model. 

Dielectric constants - The lattice contribution to the 
dielectric tensor has been calculated for both oxides in 
the fluorite, monoclinic, and Si-epitaxial structures. We 
used a standard formalism to evaluate the zero- frequency 
dielectric constant ||T^ via the frequencies of zone-center 
IR-active modes and the transverse dynamical charges. 
The vibrational modes are calculated diagonalizing the 
zone-center dynamical matrix -dF" / du^ , obtained dif- 
ferentiating by centered finite-differences (with displace- 
ments of 0.1 A) the Hellmann-Feynman force component 
a on atom i with respect to the displacement of atom 
j along direction (3. The dynamical charges are likewise 
obtained by finite-difference differentiation of the Berry- 
phase polarization with respect to atomic displace- 
ments (of typically 0.05 A). 

TABLE II: Lattice dielectric tensor for fluorite, monoclinic, 
and Si-epitaxial XO2 (the small off-diagonal elements for the 
monoclinic are not displayed for clarity), calculated using the 
dynamical charge tensor of fluorite. ef^^ is the orientational 
average measured by series capacitance in polycristalline lay- 



ers, and obtained as 




1/e"" + 1/e''" 


+ l/e 


zz 

1 




^lat 




^lat 


^avc 
= lat 


Hf02 fluorite 
IIf02 monoclinic 
Hf02 Si-epi 


27.8 
17.5 
27.6 


27.8 
15.7 
18.6 


27.8 
12.4 
24.5 


27.8 
14.9 
22.9 


Zr02 fluorite 
Zr02 monoclinic 
Zr02 Si-epi 


30.5 
24.7 
22.5 


30.5 
18.3 
71.5 


30.5 
14.6 
44.9 


30.5 
18.4 
37.0 



Since the epi-oxides were optimized without con- 
straints, they have no symmetry of practical use. The cal- 
culation of the full dynamical charge tensor for all atoms 
in the complex epitaxial (as well as the monoclinic) struc- 
ture is thus rather demanding, and currently in progress. 
In Table || we give estimates of the diagonal elements of 
the lattice dielectric tensor obtained using the dynamical 
charge tensor of the fluorite phase, which is diagonal and 
isotropic, and calculated to be .Zyf=5.20 and Zq=-2.60 
for HfOs, Z|,=5.50 and ^^=-2.75 for Zr02. Of course, 
smaller dynamical charges such as found in monoclinic 
phases ||l^, |ll| will decrease the dielectric constant, es- 
pecially the zz component. Using the monoclinic cation 
charge tensors of Refs. ||l^, |ll| and imposing the Friedel 
sum rule to obtain an average anion charge tensor, we 
estimated e^z to be 9.9 and 11.9 in monoclinic Hf02 and 
Zr02 respectively, in fair agreement with previous re- 
sults. Along with our fluorite values, also in good agree- 
ment with previous calculations, this gives us confidence 
on the reliability of our procedure. 
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FIG. 3: Orientationally averaged IR intensity spectrum 
(mode dielectric constants) of Si-epitaxial Hf02 (solid) and 
Zr02 (dashed). 



energetic standpoints. The results on the electron in- 
jection barriers are partly disappointing, as the electron 
injection barrier is much smaller that at silica/Si inter- 
faces. While insufficient for hot electrons, the barrier 
should be still acceptable for standard two-dimensional 
inversion layers, whose energy levels are at about 100 
meV above the interface triangular- well bottom pO| . 
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their surface results for Zr-0 coadsorption on Si (001) 
prior to publication. We are grateful to A. Pasquarello 
and P. Pavone for helpful suggestions. This work was 
supported in part by the European Union (INVEST 
project), and by the Parallel Supcrcomputing Initiative 
of Istituto Nazionale per la Fisica dcUa Materia. 



With reference to Table ||, for hafnia we find a re- 
duction in dielectric constant compared to fluorite both 
in the Si-epi and monoclinic phases, though the latter 
is rather more dramatic, with a more than twofold de- 
crease, in agreement with previous calculations [0. For 
zirconia, we also find a similar, approximately twofold re- 
duction of the monoclinic dielectric tensor compared to 
fluorite; notably, though, a drastic enhancement is found 
in the Si-epitaxial phase. This results from the large IR 
intensity of modes at about 90 cm^^ to 140 cm^^, as 
can be seen in Fig.^, which reports the mode dielectric 
constants ||l^ for both materials in the Si-epi phase. The 
two lower-energy modes for zirconia (dashed lines) con- 
tribute mostly to the yy component, the third to the zz 
component. The pronounced softness of Si-epi zirconia 
is presumably due to the backfolding of zone-border {X- 
point) modes. 

We carefully checked against artifacts by accurately 
reoptimizing structures and repeating phonon calcula- 
tions for different displacements. We arc confident in 
our procedure also in view of the results for the other 
phases. The single zone-center IR-active mode of fluorite 
is ti;=230 cm~^ for Hf02 and a;=258 cm^^ for Zr02; for 
the latter this agrees with recent predictions H, |l^, for 
the former the frequency is 20% lower than in Ref. 
We checked that the same results are obtained (within 
0.5% for the lattice constant and 2% for the frequency) 
with the all-electron PAW method with valence semicore 

^ . The details of the vibrational spectrum of the epi 
and monoclinic phases will be reported elsewhere, but 
we note in passing that the results for the monoclinic are 
close to previous reports . 

In conclusion, the picture of zirconia and hafnia as Si- 
gate oxides as it emerges from this work is rather en- 
couraging, certainly so from the dielectric and epitaxy- 



[1] H. R. Huff, J. Electrochem. See. 149, S35 (2002); S. 

Thompson, P. Packan, and M. Bohr, Intel Techn. J. 

Q398, p.l (2001). 
[2] D. A. Buchanan, IBM J. Res. Devel. 43, 245 (1999); G. 

D. Wilk, R. M. Wallace, and J. M. Anthony, J. Appl. 

Phys. 89, 5243 (2001). 
[3] G. E. Moore, Electronics 83, 114 (1965); SPIE 2438, 

2-17 (1995). 

[4] Y. Kim et al, IEEE Int. Electr. Dev. Meeting Techn. 

Dig. (Cat. NO.01CH37224), p. 20. 2. 1-4 (2001). 
[5] J. P. Perdew, in Electronic Structure of SoUds, P. Ziesche 

and H. Eschrig eds. (Akademie-Verlag, Berlin 1991), 

p.ll. 

[6] G. Kresse and J. Furthmiiller, Comput. Mater. Sci. 6, 15 

(1996); Phys. Rev. B 54, 11169 (1996); G. Kresse and D. 

Joubert, tbid. 59, 1758 (1999). 
[7] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 
[8] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
[9] G.-M. Rignanese, F. Detraux, X. Gonze, and A. 

Pasquarello, Phys. Rev. B 64, 134301 (2001) 
[10] X. Zhao and D. Vanderbilt, Phys. Rev. B 65, 075105 

(2002). 

[11] X. Zhao and D. Vanderbilt, Phys. Rev. B 65, 233106 
(2002). 

[12] F. Bernardini and V. Fiorentini, Phys. Rev. B 57, R9427 
(1998). 

[13] P. Blochl and C. Forst, to be published. 
[14] A. Christensen and E. A. Carter, Phys. Rev. B 62, 16968 
(2000). 

[15] M. Peressi, N. Binggeli, and A. Baldereschi, J. Phys. D 
31, 1273 (1998). 

[16] X. Zhu and S. G. Louie, Phys. Rev. B 43, 14142 (1991). 

[17] B. Kralik, E. K. Chang, and S. G. Louie, Phys. Rev. B 
57, 7027 (1998). 

[18] J. Robertson, J. Vac. Sci. Technol. B 18, 1785 (2000). 

[19] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 
1651 (1992); VASP implementation by Martijn Mars- 
man. 

[20] T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 
54, 437 (1982). 



